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Abstract. The J/-)/> is considered to be among the most important probes for the 
deconfined quark gluon plasma (QGP) created by relativistic heavy ion collisions. 
While the J/?/) is thought to dissociate in the QGP by Debye color screening, there 
are competing effects from cold nuclear matter (CNM), feed-downs from excited 
^^> ' charmonia (xc and ip') and bottom quarks, and regeneration from uncorrelated 

charm quarks. Measurements that can provide information to disentangle these 
effects are presented in this paper. 
X 

D 

^ ' 1. Introduction 

The yield of heavy quarkonia is expected to be suppressed in the QGP due to the 
^- ^ Debye screening of the color charge [T]. The J/ij^ is especially promising because of 

^ . its large production cross section and di-lepton decay channels, which make it easily 

\^ ' detected. The PHENIX experiment at RHIC is able to detect the J/ij) at midrapidity 

■"sj" . {\y\ < 0.35) via its decay to e+e~ and at forward rapidity (1.2 < \y\ < 2.2) via 

^r ' its fj,~^n~ decay. Models of J/i/' production in heavy ion collisions at RHIC energy 

contain a number of important competing effects, including modification of the J/V' 



00 

o 



X 



c^ 



in the QGP, reduced feed-down from excited charmonium states that melt just above 
the QGP transition temperature, bottom quark decay and enhancement of the yield 
due to coalescence of uncorrelated charm pairs [2]. The PHENIX Au-|-Au data at 
y/sNN = 200 GeV showed that J/tjj suppression at forward rapidity is larger than 
that at midrapidity and the suppression at midrapidity is similar to that observed by 
Vh ' NA50 at SPS in Pb-|-Pb coUisions at y/swN = 17.3 GeV [3]. However, these results are 

not well understood theoretically. Systematic study of J/ip production in heavy ion 
collisions across the entire range oiNpart is needed to disentangle the competing effects. 
PHENIX recorded Cu-|-Cu coUisions at ^/snn = 200 GeV to obtain precise data in 
the range Npart < 126, where Au-|-Au data is limited by statistics and systematic 
uncertainty and the CNM effects might be dominant [4]. Feed-downs into J/t/i in 
p + p collisions were measured at the same energy and are important in the picture of 
the sequential dissociation of quarkonia by different binding energy. Measurements of 
J/ijj in Cu-I-Cu and Au-|-Au coUisions and feed-downs into J/^ vap + p collisions are 
presented in this paper. The elliptic flow of J/i/' can set a constraint on coalescence 
of charm quarks and the measurement is reported in [S] . 
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2. Measurement of feed-downs into J/ip in p + p collisions 

An inclusive J /ip measurement in p + p collisions at ^/s=2QQ GeV has been reported 
by PHENIX [6] . Recently, PHENIX separately measured the feed-down contributions 
of Xc) i^' and bottom quarks to J /ip in p + p collisions at midrapidity {\y\ < 0.35) at 
the same energy. The preliminary results are described in this section. 

The fraction of J/tp from the Xc decay, F^^^ , was measured via the decay chain of 
Xc -^ J/'P + 7 ^- e"'"e~7 with 4.1 x 10"^ reconstructed J/ip- Branching ratios of Xcj 
to J/ip + 7 are 1.3% (xco): 36% {xci) and 20% (Xc2) Hi and the contribution of Xco 
was neglected in the analysis. The 90% confidence level upper limit obtained for _F^^ 
is 0.42. Figure [T] shows the upper limit of F^^ with results from other experiments. 
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Figure 1. The fractions of J/i/) from the Xc decay, -Fxo ™ various collision 
systems as a function of the center of mass energy, y/s. The 90% confidence 
level upper limit of -F^^. in p + p collisions at y^=200 GeV which is the PHENIX 
preliminary result is represented by the arrow. The experiment names are written 
in the figure. 



The ratio of cross sections of ^' to J/ip was measured in the e+e^ decay mode 
and the feed-down fraction of J/4> from tp' decay is 0.086 ± 0.025. The cross section 
of bottom quarks was measured by electron-hadron correlation [8] and the feed- 
down fraction of J/ip from decay of bottom and anti-bottom quarks is 0.036tQQ23 
for PT,j/ip >0 GeV/c. Branching ratios of [7] were used to obtain these feed-down 
fractions. A set of theoretically predicted feed-down fractions 9J is almost consistent 
with the PHENIX preliminary results. 



3. Measurement of J /ip in Cu+Cu and Au+Au collisions 

The final results of the J /-p measurement in Cu-l-Cu collisions are reported here and 
in [3] with the published results in Au+Au collisions ^. The numbers of observed 
J/ip in Cu-|-Cu collisions at mid and forward rapidity are about 2050 and 9000, 
respectively. Compared to the Au-|-Au data, the Cu-(-Cu data provides a more precise 
Raa measurement for Npart < 100, where the CNM effects might be dominant. 
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Figure 2. Raa vs transverse momentum pT for J/i/i production in Cu+Cu 
collisions via e+e~ decay at midrapidity (closed circles) and li~^ fJ,~ decay at 
forward rapidity (open circles). Centrality bins are written in figures. No strong 
Pt dependence of Raa is observed in Cu+Cu collisions. 
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Figure 3. a) Raa vs Npart for J/'^ production in Cu+Cu and Au+Au 
collisions at midrapidity with the prediction curves [14) with the EKS98 shadowing 
model |15| . b) The same figure at forward rapidity, c) Forward/mid rapidity 
Raa ratio, d, e, f) The same figures with the predictions curves with the nDSg 
shadowing model [16] . The breakup cross sections for the prediction curves are 
obtained from fits to the d+Au data for each rapidity bin. If the forward/mid 
rapidity Raa ratio is less than unity, it means rapidity narrowing. 
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Figure [H shows Raa as a function of px for each centrality bin in Cu+Cu coUisions 
at mid and forward rapidity. Suppression by a factor of two is seen at both mid and 
forward rapidity in the most central coUisions. No strong px dependence of Raa is 
observed. Analysis is ongoing to extend Raa beyond pt—5 GeV/c, where the STAR 
experiment has a moderate acceptance [10] . In central Au+ Au collisions, Raa slightly 
increases with px at forward rapidity and Raa is flat within errors below pT=5 GeV/c 
at midrapidity. 

Figure [3] shows Raa at mid and forward rapidity and the forward/mid rapidity 
Raa ratio in Cu+Cu and Au+Au collisions. The Cu+Cu and Au+Au data agree 
well in the overlap region. The results of the J/ip measurement in d+Au collisions, 
which is the primary tool to investigate the CNM effects, were updated and are 
reported in [TTl [1^. PHENIX has interpreted the d+Au data and extrapolated 
prediction of the CNM effects to heavy ion collisions based on calculations of R. Vogt 
including nuclear shadowing and nuclear breakup [T31 HI] • While the nuclear breakup 
cross sections of J/i^, crbreakup, does not vary with rapidity in the framework of the 
calculations, PHENIX treated (Jf^reakup as a rapidity dependent "effective" parameter 
and determined it by fitting to the d+Au data [4]. The theoretical prediction curves 
with ±lcr bands obtained for the rapidity dependent "effective" (Tbreakup with two 
nuclear shadowing models (EKS98 [15] and uDSg [16]) are also shown in figure [Sj 
Since the obtained Ubreakup at forward rapidity is larger than one at midrapidity, 
the expected CNM effects from the PHENIX "effective" model at forward rapidity 
are larger than those at midrapidity and this leads to rapidity narrowing. The J/ip 
suppression beyond CNM effects seems to start at Npart ~ 200 (100) at mid (forward) 
rapidity. As shown in figurejSjc) and f), rapidity narrowing in central Au+Au collisions 
is consistent with the expected CNM effects. 

While the current uncertainty of the CNM effects is large, it is expected to be 
reduced by the d+Au collision data collected in 2007-2008 whose statistics are 30 
times larger than the currently available d+Au statistics. 
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